Sudden transition between classical and quantum decoherence in dissipative cavity 

QED and stationary quantum discord 
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We investigate the phenomenon of sudden transition between classical and quantum decoherence 
in the study of quantum discord for a dissipative cavity QED system, which consists of two nonin- 
teracting two-level atoms, each trapped in a dissipative cavity. It is found that the quantum discord 
between two atoms, which are prepared initially in the X-type quantum states, is not destroyed 
by the dissipation of the cavities for a finite time interval and the stationary quantum discord can 
arise in the interaction of atoms with cavities as the time approach to infinite. The transition time 
is sensitive to the initial state parameter of the two atoms and the mean photon number of the 
coherent field. Interestingly, the quantum discord between the two atoms is completely unaffected 
by the dissipation of the cavities if we choose the suitable value of the ratio, which depends on the 
decay rate of two cavities and the atom-field coupling constant. 



I. INTRODUCTION 

Quantum correlations play an important role in quan- 
tum information and computation and have become a 
hot topic of intense research in recent years. Entangle- 
ment is a special kind of quantum correlation and has 
been recognized as an essential resource for many opera- 
tions in quantum information processing [1-4]. However, 
the entanglement is not the only type of quantum cor- 
relation and there exist quantum tasks that display the 
quantum advantage without entanglement [5-7]. It has 
been demonstrated both theoretically [8-10] and experi- 
mentally [11] that other nonclassical correlation, namely, 
quantum discord [12] can be responsible for the computa- 
tional speedup for certain quantum tasks. Quantum dis- 
cord, introduced in Ref. [12], is defined as the difference 
between the quantum mutual information and the clas- 
sical correlation and is nonzero even for separate mixed 
states. Therefore, the quantum discord may be regarded 
as a more general and fundamental resource for quantum 
information processing. 

On the other hand, a real quantum systems will un- 
avoidably be influenced by surrounding environments. 
The interaction between the quantum system and its en- 
vironment leads to a rapid destruction of quantum co- 
herence, which is the main problem for the realization 
of quantum information processing. Therefore, it is very 
important to understand the dynamics of quantum cor- 
relations of open systems and find potential applications 
in quantum information theory. Many efforts have been 
devoted to the study of the dynamics of quantum cor- 
relation under various decoherence channels [13-16]. It 
has been shown that the discord is more robust than en- 



xujb@zju.edu.cn 



tanglement under the Markovian environments [14] and 
vanishes only at some time points under non-Markovian 
environment [16]. 

Recently, the dynamics of quantum and classical cor- 
relations have been studied in the presence of nondissipa- 
tive decoherence and the phenomenon of the sudden tran- 
sition from classical to quantum decoherence in a finite 
time interval has been reported [17]. Some of these phe- 
nomena have been observed in the recent experiment [18] . 
In this paper, we investigate the phenomenon of sudden 
transition between classical and quantum decoherence in 
the study of quantum discord for the dissipative cavity 
QED system, which consists of two noninteracting two- 
level atoms, each trapped in a dissipative cavity. Firstly, 
we assume that the two atoms are prepared initially in 
the X-type quantum states and show that the quantum 
discord of two atoms is not destroyed by the dissipation 
of the cavities for a finite time interval and revivals to a 
stable value after damping oscillation, which means that 
there is a stationary quantum discord between the two 
atoms as the time approach to infinite. This implies that 
the initial quantum discord of two atoms can be partially 
preserved even when they are put into the two spatially 
separated dissipative cavities, respectively. Furthermore, 
we can see clearly that the transition time depends on the 
parameter of the initial states of the two atoms and the 
mean photon number of the coherent field. Particularly, 
it is interesting to point out that the quantum discord of 
two atoms is completely unaffected by the dissipation of 
the cavities if we choose the suitable ratio, which depends 
on the decay rate of two cavities and the atoms-field cou- 
pling constant. Then, we find that the sudden transition 
phenomenon does not appear in this system if the two 
atoms are initially in the Werner state. Instead, the sta- 
tionary quantum discord of the two atoms still exists in 
the long-time regime. It is worth noting that the amount 
of the stationary quantum discord between two atoms 
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FIG. 1. This is the schematic diagram of the system which 
is investigated in the present paper. There is no interaction 
between two atoms and no communication between two cav- 
ities. 



can be enhanced by increasing the value of the ratio. 



be read as 



dp(t) 
dt 



-i[H eS ,p{t)]+Dp(t) 

-i[H eS , p(t)} + D AP (t) + D B p(t), (3) 



where the superoperators Da and Db represent the losses 
in the cavity A and B. At zero temperature, the super- 
operators Da and Db can be written as 



Da = k(2aA ■ a\ - a' A a A a A a A ), 

Db = k(2a B ■ a B - a B a B a B a B ), 



(4) 



where k is the cavity decay rate. For simplicity, we as- 
sume the decay rates of two cavities are equal. In the 
interaction picture, the master equation takes the form 



dp(t) 
dt 



-i[V,p(t)]+D A p(t) + D B p(t), (5) 



II. THE DYNAMICS OF TWO ATOMS IN 
DISSIPATIVE CAVITIES 



We consider the system consisting of two noninteract- 
ing two-level atoms, each trapped inside a dissipative 
cavity (see Fig. 1). The Hamiltonian describing the in- 
teraction between atoms and cavities can be written as 
(h=l) 

H =^{cr A +a B )+u{a A a A +b\b A )+g ]T (at CT r +a l a+) 7 

i=A.B 

(1) 

where g is the atom-field coupling constant, a 1 and a 
are the creation and annihilation operator of the single- 
mode cavity field, and a z is the atomic inversion oper- 
ator, o~l = \e)i(g\ (a," = |<7)»(e|) is the atomic spin flip 
operators. The symbols |e) and \g) refer to the excited 
and ground states of the two-level atom. cj and w are 
the atomic transition frequency and cavity frequency, re- 
spectively. It is clear that there is no interaction between 
subsystem A and subsystem B, which means that the 
evolution of each subsystem is independent. In the dis- 
persive approximation, the Hamiltonian of Eq. (1) can 
be rewritten as [19] 

#eff = y(ol + °b) + u(a\a A + b\b A ) + fl[(a A a A + 1) 

\e) A (e\ - a A a A \g) A (g\] + tt[(a B a B + 1) 
\e) B (e\- a B a B \g) B (g\], (2) 

with A = ojq — oj and fl = 5 2 /A. 

Next, we investigate the dynamical evolution of two 
noninteracting two-level atoms interacting with two dis- 
sipative cavities by making use of the master equation. 
In the dispersive approximation, the master equation can 



with 

v = n [K tfl « + l )\ e )M - «\<*\9)M- (6) 

i=A,B 

This equation can be solved by using the superoperator 
method [19-21]. We assume that the two cavities are pre- 
pared initially in the coherent states \a\) and {012), and 
the two atoms are prepared in a class of state with the 
maximally mixed marginals, which can be described by 
the three-parameter X-type density matrix in the basis 
{\ee),\eg),\ge),\gg)} as 



1 

Patom(O) = -(I + C i°A ® <*b) 

I 1 + c 3 ci - c 2 

1 1 - c 3 ci + c 2 

~ 4 ci + c 2 1 - c 3 

V ci - c 2 1 + c 3 



(7) 



where Cj (0 < |cj| < 1) are the real numbers , / is the 
identity operator of the total system, <j\ and a B are the 
Pauli operators of atom A and atom B, respectively. For 
simplicity, we set a\ = a 2 = a. Then, the density matrix 
of the whole system at t — is 



/ l + c 3 



P(0) = s 











ci - c 2 \ 



1 — c 3 ci + c 2 
ci + c 2 1 — c 3 

V ci - c 2 1 + c 3 J 

ol)a(ol\ ® |a) B (a|. (8) 



It is not difficult to find that the density matrix of the 
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system at time t is 

1 + C 3 

p(t) = — - — \ee) AB (ee\\a + {t)a + (t)) AB (a + (t)a + {t)\ 

1 - c 3 

H —\eg) AB(eg\\a + (t)a-(t)) AB (a + (t)a-(t)\ + 

1 — C3 1 ~|- C3 

—4— \ge) A B(ge\\a-(t)a + (t)) A B(a-(t)a + {t)\ + 

\gg) A B(gg\\a-(t)a_(t)) AB (a_(t)a_(t)\+{^^f(t) 2 

\ee) AB (gg\\a + (t)a + (t)) AB (a-(t)a-(t)\ + C ^±^\f(t)\ 2 
\eg) AB (ge\\a + (t)a_(t)) AB (a_(t)a + (t)\ + H.C.}, (9) 
with 

f(t) = exp{-ifit + \a\ 2 (e- 2kt - 1)} 

|a±(i)) = \aexp-(k±in)t). (10) 
where if.C represents the Hermitian conjugate. 

III. THE DYNAMICS OF QUANTUM DISCORD 
OF TWO ATOMS 

In this section, we investigate the dynamics of quan- 
tum discord of two atoms which are trapped in the two 
spatially separated and dissipative cavities, respectively. 
For a two-qubit quantum system, the quantum discord 
qualifying a measure of quantum correlation, is defined as 
the difference between the quantum mutual information 
and the classical correlation, 

Q( PAB ) = T{p AB ) - C(p AB ), (11) 

where I{p AB ) is the total correlation of two subsystem 
and can be expressed as 

l(pAB) = S(p A ) + S(p B )-S(p AB ), (12) 

Here, S(p) = —tr(p\og 2 p) is the von Neumann en- 
tropy, p A and p B are the reduced density matrices of 
p AB . Besides, C{p AB ) is the classical correlation between 
two subsystem A and B, which is defined as the maxi- 
mum information one can obtain form A by performing 
a perfect measurement on B. As discussed in Ref. [22], 
the classical correlation is described as 

C( PAB ) = mzx{S( PA ) - S{p AB \{B k }), (13) 

{B k } 

where {B k } is a complete set of projectors preformed 
locally on subsystem B. S(p A B\{B k }) = J2 k PkS{p k ) 
is the based-on-measurement quantum conditional en- 
tropy, (p AB \{B k }) := p k = l/pk(I ® B k )p AB (I <g> B k ) 
is the conditional density operator and p k = tr( AB ) [{I ® 
B k )p AB (I <8> B k )] is the probability. 

Next, we begin to study the quantum discord dynamic 
properties of two atoms. Tracing over the degrees of the 



freedom of cavity fields in the Eq. (9), we can obtain the 
reduced density matrix of atom A and atom B, 

PA B (t) = ^-^-\ee)(ee\ + ^—^-\eg){eg\ + 

\ge)(ge\ + ±^\gg){gg\ + { C -l^l { f(t) X {t)f 

\ee)(gg\ + C -l±^\f(t) X (t)\ 2 \eg)(ge\ + H.C.}, (14) 

where x(i) = (ct-(t)\a + (t)} and H.C. is the Hermitian 
conjugate. 

The eigenvalues of the reduced density matrix p AB (t) 
in Eq. (14) can be derived as, 

Ai,2 = \[(l-c 3 )±\f(t) X (t)\ 2 (c 1+ c 2 )}, 

A 3 ,4 = \[(l + c 3 )±\f(t) X (t)\ 2 ( Cl -c 2 )]. (15) 

It is not difficult to find from Eq. (14) that the reduced 
density matrix of each subsystem is maximally mixed, 
which means that p A (t) = p B {t) = 1/2. Consequently, 
the von Neumann entropy S(p A (t)) = S(p B (tj) = 1. 
Then the quantum mutual information between two sub- 
systems is 

4 

I{pA B (t)) = 2 + ^A i log 2 A i . (16) 

i=l 

In order to calculate the classical correlation C(p AB ), 
we choose the complete set of projectors {B k = \9 k ){0 k \, 
(k = 1,2)} to measure the subsystem B, where the two 
orthogonal projectors are defined by, 

|0i) =cos0|e) + e^sin%), 

|02> =e"^ single) -cos%), (17) 

with the parameters 9 and <f) vary from to 2ir, respec- 
tively. After the measurement {B k , (k = 1, 2)}, the prob- 
ability pi = P2 = 1/2 and the reduced matrices of sub- 
system A can be obtained as follows, 

p\{t) = i[(l + C3 )cos 2 #+(l- C3 )sin 2 6}\e) A (e\ 

+ i[(l - c 3 ) cos 2 9 + (1 + c 3 ) sin 2 9} \g) A (g\ + 

\{[(ci - c 2 )e^f 2 (t) X 2 (t) cos 9 sin 9 + ( Cl + c 2 ) 
e-^\f(t) X (t)\ 2 cos9sm9]\e) A (g\ + H.C.}, 
p\{t) = i[(l-c 3 )cos 2 #+(l + c 3 )sin 2 9}\e) A (e\ 

+ i[(l + c 3 ) cos 2 9 + (1 - c 3 ) sin 2 9} \g) A (g\ + 

\{[-{ci - c 2 y*f 2 (t) X 2 (t) cos 9 sin 9- (ci + c 2 ) 
e~ t,t '\f(t)x(t)\ 2 cos9sm9}\e) A (g\ +H.C.}, (18) 
where H.C. denotes the Hermitian conjugate. 
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FIG. 2. Time evolution of quantum mutual information (dot- 
ted line), classical correlation (dashed line) and quantum dis- 
cord (solid line) of two atoms as a function of the dimension- 
less scaled time Sit with k/Sl = 0.05 and C3 = 0.6. Vertical 
dashed line corresponds to the transition time Sit. Panel (a): 
a = 1.2. Panel (b): a = 0.8. 



The eigenvalues of the above reduced density matrix 
Pa V s = 1>2) can be calculated as (note that the eigen- 
values of two reduced density matrices are equal), 



.(*) _ 



1 



- 2< 1:: V), 



(19) 



where the parameter r\ is defined by 



1 



with 



V = {ci cos' 26 +-\f(t) X (t)\ i [2(c 2 1 +c 2 2 ) 
+2{c\ - c|) sin(2</> + if)] sin 2 2<9} 1/2 , 

/W(*) + /*V 2 



(20) 



sm ip 
cos ip 



2|/(%(*)| 2 ' 

»(/WW - rv 2 ) 

2|/(%(*)| 2 



(21) 



Consequently, the von Neumann entropies S(p A (t)) 
and S(p A (t)) are given by 

S(p A (t)) = S{ P \{t)) 



1 — 77 , 1 — 77 1 + 77, 1 + 7/ 
-1T~ l0 §2 a !0g 2 "Tr - 
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FIG. 3. The quantum discord of two atoms is plotted as a 
function of the dimensionless scaled time Sit and parameter 
C3 with a — 1 and k/Sl = 0.05. 



Using Eq. (13), the classical correlation can be written 



as 



C(p AB {t)) = max{5(pA) - SVabK-BaJ) 

= 1 -^(y\pkS{pk)} 

&,<p — 
k 

= 1 — imn[R(rf)}. 

B,(f> 



(23) 



Since the function R(i]) is a monotonically decreasing 
function, we can obtain the minimal value of R{rf) by 
choosing the suitable parameters 9 and <fi to ensure that 
the parameter 77 defined in Eq. (18) is maximal. From 
the Eq. (18), it is not difficult to find that there exists an 
inequality 



77 < {c 2 cos 2 29 



\f(t)x(t)\ 4 



[2(c 



-2|c 2 ~c 2 |]sin 2 20} 1 / 2 

| C3 | ii\c 3 \>W(t) 

< W(t) if leal < W(t) 



with 



(24) 



(25) 



Combining Eqs. (23) and (24), we can rewrite the clas- 
sical correlation as 



C{p AB {t)) = 2^ 5 



(22) where m(i) = max{|c 3 |, W(*)}. 



logati + c-ijM*)]- 

(26) 
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FIG. 4. Time evolution of quantum mutual information (dot- 
ted line), classical correlation (dashed line) and quantum dis- 
cord (solid line) of two atoms as a function of the dimension- 
less scaled time Qt with k/Q = 0.05, a = 0.8 and C3 = 0.6. 
Panel (a): < Q < 10. Panel (b): < Q < 50. 



FIG. 5. Time evolution of quantum mutual information (dot- 
ted line), classical correlation (dashed line) and quantum dis- 
cord (solid line) of two atoms as a function of the dimension- 
less scaled time Qt. Panel (a): a = 0.8, c a = 0.7 and k/Q = 2. 
Panel (b): a = 1.2, c 3 = 0.6 and k/Q = 3. 



Then the quantum discord between two atoms can be 
expressed as 

Q(pAB(t)) = T{ PAB {t)) - C( P AB{t)) 
4 

= 2 + Ai log 2 A, 

»=i 

+ m(t) log 2 [l + (-l)M*)h (27) 

We first focus on the X-type quantum states with state 
parameters c\ = 1, c 2 = — c 3 and | C3 1 < 1. Substituting 
these initial conditions into the Eqs. (16), (26) and (27), 
the explicit expressions of the quantum mutual informa- 
tion, classical correlation and quantum discord can be 
obtained as, 

APA B {t)) = \{l + c 3 ) log 2 (l + c 3 ) + ^(1 - c 3 ) 

log 2 (l - c 3 ) + \{l + \f(t) X (t)\ 2 ) log 2 (l + \f(t) X (t)\ 2 ) 

+ ^(l-|/(%(i)| 2 )log 2 (l-|/(%(i)| 2 ), 

C( PAB (t)) = E 1 + m{t) log 2 [l + (-l)M*)]. 

i=i 

Q{pA B it)) - l(p AB {t)) - C( PAB {t)). (28) 
where m(£) = max{|c 3 |, \ f{t)x(t)\ 2 }- 



In Fig. 2, we plot the quantum discord (solid line), 
the classical correlation (dashed line) and the quantum 
mutual information (dotted line) as a function of the di- 
mensionless scaled time ilt for two different values of a 
with k/Q. = 0.05 and c 3 = 0.6. It is shown that the 
phenomenon of the sudden transition between classical 
and quantum decoherence appears in this system within 
a short interaction time. Comparing the panel (a) with 
panel (b), we find that the quantum discord between two 
atoms is not destroyed by the dissipation of the cavities 
for a finite time interval and the transition time fit can 
be tuned by adjusting the value of the parameter a. 

In order to illustrate the dependence of transition time 
on the initial states, the quantum discord Q(p) is dis- 
played as a function of the dimcnsionless scaled time fit 
and the parameter c 3 with a = 1 and k/Q = 0.05 in 
Fig. 3. We can see that the quantum discord is a constant 
over the time interval Qt < Qt and can be increased by 
decreasing the value of c 3 . This consequence is in agree- 
ment with the conclusion which is obtained by Ref. [17] . 

Since the decay of quantum and classical correlations 
in this system corresponds to hyperexponential decay, 
the long-time behavior is different from the situation of 
exponential decay of Ref. [17]. In Fig. 4, we plot the 
correlations as a function of the dimensionless scaled time 
Qt (0<Qt< 50) with a = 0.8, k/Q = 0.05 and c 3 = 0.6. 
It is quite clear that the quantum correlation revivals 
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FIG. 6. Time evolution of quantum discord of two atoms as 
a function of the dimensionless scaled time fit with a = 0.5 
and k/fl = 0.05 for r = 0.5 (dotted line) and r = 0.9 (solid 
line). Panel (a): < fit < 10. Panel (b): < fit < 50. 



FIG. 7. Time evolution of quantum discord of two atoms as 
a function of the dimensionless scaled time fit with r — 0.7 
and k/fl — 0.05 for a — 0.8 (dotted line) and a — 0.5 (solid 
line). Panel (a): < fit < 10. Panel (b): < fit < 50. 



to a stable value after damping oscillation, which means 
that the initial quantum discord of two atoms can be 
partially preserved even when they are put into the two 
spatially separated dissipative cavities, respectively. The 
amount of stationary quantum discord depends on the 
initial state of the two atoms and the parameter a. 

The influence of the ratio k/fl on the dynamics evo- 
lution of the quantum discord (solid line), classical cor- 
relation (dashed line) and quantum mutual information 
(dotted line) is displayed in Fig. 5. Comparing the Fig. 5 
with the Fig. 4, we can see that the quantum correlation 
of two atoms is completely unaffected by the decoher- 
ence of the cavities if we choose the suitable ratio k/fl. 
The reason is that the cavity fields deplete quicker with 
increasing the dissipation parameter k, the chance for 
atoms to interact with the cavity fields decreases. This 
result may have some applications in the quantum infor- 
mation processing and quantum memory. 

In the following, we consider the class of initial states 
for which c\ = c 2 = C3 = —r, with r G [0,1]. In this 
situation, the initial state p a tom(0) turns out to be the 
Werner state [23] 

Patom(0) = (l-r)|+r|*-)<*-|, (29) 



with 



[*-) = -j={\eg) (ge\). 



(30) 



Using the similar procedure, we can obtain the quan- 
tum mutual information as 



ApABit)) = 2 + ^A 2 log 2 Ai, 



i=i 
l-r 



Ai — A2 

r 

A 3 = A 4 = i(l + r±2r|/(t) X (t)| 2 ), (31) 



the classical correlation is given by 



2 

C( P AB(t)) = £ 1 + ( ~ 1)M * } log 2 [l + (-l) J 'n(t)], 



n(t) = max{r,r|/(t)xW| 2 }, 
and the quantum discord is 

Q(pA B {t)) = I(pA B (t)) - C(p AB (t)) 



(32) 



log 2 [l + r+(-l)^2r|/(t) X (t)| 2 ]-^ 



1 + (-l) J 'n(t) 



log 2 [l + (-l)''n(i)]. 



(33) 



In Fig. 6, we plot the quantum discord Q(p) as a func- 
tion of the dimensionless scaled time fit with a = 0.5 
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FIG. 8. Time evolution of quantum discord of two atoms as 
a function of the dimensionless scaled time fii for k/Q, = 0.1 
(dotted line) and k/Q = 1 (solid line). Panel (a): r = 0.7 and 
a = 0.8. Panel (b): r = 0.9 and a = 0.5. 



and k/Q = 0.05 for r = 0.5 (dotted line) and r = 0.9 
(solid line). It is shown that the phenomenon of sudden 
transition does not appear in this situation. Instead, the 
quantum discord revivals to a stable value after damping 
oscillation, which means that there is a stationary quan- 
tum discord between two atoms as the time approach to 
infinite. Furthermore, we also find that stationary quan- 
tum discord of two atoms increases with increaseing the 
parameter r. 

The evolution of the quantum discord of two atoms is 
plotted with r = 0.7 and k/Q = 0.05 for a = 0.8 (dotted 
line) and a = 0.5 (solid line) in Fig. 7. It is obvious 
that there is a stationary quantum discord between two 
atoms and can be increased by decreasing the value of 
the parameter a. 

In order to show the influence of ratio k/Q on the 
behavior of quantum discord of two atoms, we plot the 
quantum discord as a function of the dimensionless scaled 
time Qt for two different values of the ratio k/Q in Fig. 8. 



We can see clearly from Fig. 8 that for the different initial 
conditions, the quantum discord of the stationary state 
can be enhanced by increasing the value of the ratio k/Q. 



IV. CONCLUSIONS 

In this paper, we investigate the phenomenon of sud- 
den transition between classical and quantum decoher- 
ence in the study of quantum discord for the dissipativc 
cavity QED system, which consists of two noninteract- 
ing two-level atoms, each trapped in a dissipative cavity. 
Firstly, we assume that the two atoms are prepared ini- 
tially in the X-type quantum states with state parame- 
ters ci = 1, c 2 — — C3 and | C3 j < 1. It is shown that the 
quantum discord of two atoms is not destroyed by the 
dissipation of the cavities for a finite time interval and 
partially recovers its initial values for a long interaction 
time, which means that there is a stationary quantum 
discord between the two atoms as the time approach to 
infinite. This demonstrates that the initial quantum dis- 
cord of two atoms can be partially preserved even when 
they are put into the two spatially separated dissipative 
cavities, respectively. Furthermore, we notice that the 
transition time Qt depends on the parameter of initial 
states of the two atoms and the mean photon number 
of the coherent field. Particularly, it is interesting to 
point out that the quantum discord of two atoms are 
completely unaffected by the dissipation of the cavities 
if we choose the suitable value of the ratio k/Q, which 
depends on the decay rate of two cavities and the atoms- 
field coupling constant. Then, we find that the sudden 
transition phenomenon does not appear in this system if 
the two atoms are initially in the Werner state. Instead, 
the stationary quantum discord of the two atoms still ex- 
ists in the long-time regime. It is worth noting that the 
amount of the stationary quantum discord between the 
two atoms can be enhanced by increasing k/Q. These 
results may have potential applications in the quantum 
information processing and quantum memory. 
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